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Liquids in contact with solids are submitted to intermolecular forces making 
liquids heterogeneous and stress tensors are not any more spherical as in homo- 
geneous bulks. The aim of this article is to show that a square-gradient func- 
tional representing liquid-vapor interface free energy corrected with a liquid 
density functional at solid surfaces is a well adapted model to study structures 
of very thin nanofilms near solid walls. This result makes it possible to study 
the motions of liquids in nanolayers and to generalize the approximation of 
lubrication in long wave hypothesis. 

Keywords: Nanolayers, disjoining pressure, thin flows, approximation of lubri- 
cation. 



1. Introduction 

At the end of the nineteenth century, the fluid inhomogeneity in liquid- 
vapor interfaces was taken into account by considering a volume energy de- 
pending on space density derivative."'^ This van der Waals square-gradient 
functional is unable to model repulsive force contributions and misses 
the dominant damped oscillatory packing structure of liquid interlayers 
near a substrate wall.^ Furthermore, the decay lengths are correct only 
close to the liquid-vapor critical point where the damped oscillatory struc- 
ture is subdominant.'^ In mean field theory, weighted density-functional 
has been used to explicitly demonstrate the dominance of this structural 
contribution in van der Waals thin films and to take into account long- 
wavelength capillary-wave fluctuations as in papers that renormalize the 
square-gradient functional to include capillary wave fluctuations.'* In con- 
trast, fluctuations strongly damp oscillatory structure and it is mainly for 
this reason that van der Waals' original prediction of a hyperbolic tangent is 
so close to simulations and experiments.^ The recent development of exper- 
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imcntal technics allows us to observe physical phenomena at length scales 
of a few nanometers.^ To get an analytic expression in density- functional 
theory for liquid film of a few nanometer thickness near a solid wall, we 
add a liquid density-functional at the solid surface to the square-gradient 
functional representing closc^ly liquid- vapor interface free energy. This kind 
of functional is well-known in the literature.^ It was used by Cahn in a phe- 
nomenological form, in a well-known paper studying wetting near a critical 
point.* An asymptotic expression is obtained in^ with an approximation 
of hard sphere molecules and London potentials for liquid-liquid and solid- 
liquid interactions: we took into account the power-law behavior which is 
dominant in a thin liquid film in contact with a solid. 

For fluids submitted to this density-functional, we recall the equation of 
motion and boundary conditions. We point out the definition of disjoining 
pressure and analyze the consequences of the model. Finally, we study the 
motions in liquid nanolayers; these motions are always object of many de- 
bates. Within lubrication and long wave approximations, a relation between 
disjoining pressure, viscosity of liquid, nanolayer thickness variations along 
the layer and tangential velocity of the liquid is deduced. 

2. The density- functional 

The free energy density-functional of an inhomogeneous fluid in a domain 
O of boundary dO is taken in the form 

F = [ [ [ sdv + [ [ ipds. (1) 

J J Jo J JdO 

The first integral is associated with square-gradient approximation when 
we introduce a specific free energy of the fluid at a given temperature 0, 
e = s{p, (3) as a function of density p and /3 = (grad p)^. Specific free energy 
s characterizes together fluid properties of compressibility and molecular 
capillarity of liquid- vapor interfaces. In accordance with gas kinetic theory, 
A = 2pe'ij{p,(3) is assumed to be constant at given temperature^^ and 

pe = pa{p) + ^{gv&d pf, (2) 

where term (A/2) (grad p)^ is added to the volume free energy pa{p) of a 
compressible fluid. Speciflc free energy a enables to connect continuously 
liquid and vapor bulks and pressure P(p) = p'^a'p{p) is similar to van dcr 
Waals one. Near a solid wall, London potentials of liquid-liquid and liquid- 
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solid interactions are 

ipil = — ^ , when r > ai and ipii = oo when r < ai , 

Vis = ——^ , when r > S and ipis = oo when r < 6 , 

where cu and c/s are two positive constants associated with Hamakcr con- 
stants, CT; and as denote fluid and solid molecular diameters, 6 = i( ai + 
as) is the minimal distance between centers of fluid and solid molecules. 
Forces between liquid and solid have short range and can be described 
simply by adding a special energy at the surface. This is not the entire 
interfacial energy; another contribution comes from the distortions in the 
density profile near the wall.^'^^ For a plane solid wall (at a molecular scale), 
this surface free energy is 

^^(P) = -71P+ ^72P^- (3) 
Here p denotes the fluid density value at the wall; constants 71, 72 are 
positive and given by relations 71 = ^2S^mim ^^"''^ ^ 126^171^ ' ^^^^^ 
mi and nis denote respectively masses of fluid and solid molecules, psoi is 
the solid density.^ Moreover, we have A 



3(7; mf 

We consider a horizontal plane liquid intcrlaycr contiguous to its vapor 
bulk and in contact with a plane solid wall (S); the z-axis is perpendicular 
to the solid surface. The liquid film thickness is denoted by h. Conditions 
in vapor bulk yield gradp = and Ap = 0. Another way to take into 
account the vapor bulk contiguous to the liquid interlayer is to compute 
a density- functional of the complete liquid-vapor interlayer by adding a 
supplementary surface energy t/j on a geometrical surface (E) a,t z = h to 
volume energy (2) in liquid intcrlaycr (L) and surface energy (3) on solid 
wall (>S').^'^ This assumption corresponds to a liquid interlayer included 
between z = and z = h, a. liquid- vapor interface of a few Angstrom 
thickness assimilated to surface z = h and a vapor layer included between 
z = h and z = 00. Due to small vapor density, let us denote by tp the 
surface free energy of a liquid in contact with a vacuum, 

where 74 ~ 72 and p is the liquid density in a convenient point inside the 
liquid-vapor interface. Density-functional (1) of the liquid-vapor layer gets 
the final form 

F= s dv+ (f>ds+ ^ ds 

J J J(L) J J{S) J Ji^) 
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3. Equation of motion and boundary conditions 

In case of equilibrium, functional F is minimal and yields the equation of 
equilibrium and boundary conditions. In case of motions we simply add the 
inertial forces pF and the dissipative stresses to the results. ^^"^^ 

3.1. Equation of motion 

The equation of motion is^*^'^^ 

p r = div (ct + CTv) — p grad O , (5) 

where T is the acceleration, fl the body force potential and a the stress 
tensor generalization 

(T = — p 1 — A grad p eg) grad p, 

with p — p^e'p — p div (A grad p). The viscous stress tensor is = 
Ki ( tr I?) 1 + 2 K2 D where D denotes the velocity strain tensor; m and 
K2 are the coefficients of viscosity. 

For a horizontal layer, in an orthogonal system of coordinates such that the 
third coordinate is the vertical direction, the stress tensor a of the thin film 
takes the form : 

„ A [dpV ^ d'^p 

A / dp\ ^ d'^p 
^^ = -^-2UJ ^''d^ 
Let us consider a thin film of liquid at equilibrium {gravity forces are ne- 
glected but the variable of position is the ascendant vertical). The equation 
of equilibrium is : 

div cr = (6) 
Eq. (6) yields a constant value for the eigenvalue a^, 

where denotes pressure P{pvb) in the vapor bulk of density p^^ bounding 
the liquid layer. Eigenvalues ai,a2 are not constant but depend on the 
distance z to the solid wall.^^ At equilibrium, Eq. (5) yields:^'' 

grad [p(p)-AAp]=0, (7) 

where p is the chemical potential at temperature 9 defined to an unknown 
additive constant. The chemical potential is a function of P (and 0) but it 
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can be also expressed as a function of p (and 9). We choose as reference 
chemical potential jio = l^o{.p) null for bulks of densities pi and pv of phase 
equilibrium. Due to Maxwell rule, the volume free energy associated with 
Po is ffo(p) — Po where Po = P(,pi) = P{pv) is the bulk pressure and go{p) = 
Jp A^o (p) dp is null for liquid and vapor bulks of phase equilibrium. The 
pressure P is 

Pip) = ppo{p)-9o{p) +Po- (8) 

Thanks to Eq. (7), we obtain in the fluid and not only in the fluid interlay er, 

Po{p) - AAp = Po{pb), 

where Po{pb) is the chemical potential value of a liquid mother bulk of 
density pi, such that Po{pb) = l^oiPvb)^ where py^ is the density of the vapor 
mother bulk bounding the layer. We must emphasis that P{pb) and P{pvb) 
are unequal as for drop or bubble bulk pressures. Density pj, is not a fluid 
density in the interlayer but density in the liquid bulk from which the 
interface layer can extend (this is the reason why Derjaguin used the term 
mother liquid,^"^ page 32). In the interlayer 

d'^p 

>^-p2 = Pb{p), with pi,{p) = po{p) - Po{Pb) (9) 



3.2. Boundary conditions 

Condition at the solid wall (5) is associated with Eq. (3)^^ : 

where n is the external normal direction to the fluid; Eq. (10) yields 

Condition at the liquid- vapor interface (S) is associated with Eq. (4): 

I), <") 

' \z = h 

Eq. (11) defines the film thickness by introducing a reference point inside 
the liquid- vapor interface bordering the liquid interlayer with a convenient 
density bX, z = h}^ 

We must also add the classical surface conditions on the stress vector asso- 
ciated with the total stress tensor a + to these conditions on density. 
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4. The disjoining pressure for horizontal liquid films 

We consider fluids and solids at a given tem,perature 6. The hydrostatic 
pressure in a thin liquid interlayer included between a solid wall and a vapor 
bulk differs from the pressure in the contiguous liquid phase. At equilibrium, 
the additional pressure interlayer is called the disjoining pressure}^ The 
measure of a disjoining pressure is either the additional pressure on the 
surface or the drop in the pressure within the mother bulks that produce 
the interlayer. The disjoining pressure is equal to the difference between the 
pressure P„j_ on the intcrfacial surface (pressure of the vapor mother bulk 
of density p^^) and the pressure P}, in the liquid mother bulk (density pi,) 
from which the interlayer extends : 

n(/i) - p,, - Pfc. 

li gh{p) = 9o{p) - 9o{Pb) - IJ'o{Pb){p- Pb) denotes the primitive of /Xb(p) nuU 
for pb, we get from Eq. (8) 

n(Pfe) = -9biPv,), (12) 

and an integration of Eq. (9) yields 

The reference chemical potential linearized near pi {respectively py) is 

Po{p) = —(p — pi) {respectively p.o{p) = —{p — pv) ) where c; {respectively 

Pi Pv 
Cy) is the isothermal sound velocity in liquid bulk pi {respectively vapor hulk 

Pv) at temperature 9}^ In the liquid and vapor parts of the liquid-vapor 

film, Eq. (9) yields 

dp p (? p 
\-r^ = —{p-pb) (liquid) and \-— = -^{p-p) (vapor). 
az^ pi az^ pv 

The values of p.o{p) are equal for the mother densities pv^ and ph, 

(? 

— {Pb-Pi)= Po{Pb) = Po{Pvh) = —{Pvk - Pv), and consequently, 

Pi Pv 

(. , 4 {Pb - Pi) 

Pv^= Pv\^ + -\ 

\ cl Pi 

In liquid and vapor parts of the liquid-vapor interlayer we have, 

9o{p) ^-^{p- Pif (hquid) and go{p) = ^{p - Pvf (vapor). 

Ipi Ipv 
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From definition of gb{p) and Eq. (12) we deduce the disjoining pressure 



nipt,) = -^{pi - Pb) 
2pi 



Pi 



(14) 



Due to pj 2 + ^ < p; + Pb, we get Il{pb) « ^{pf - pi). 



Pi J "^Pi 
Now, wc consider a film of thickness h: the density profile in the liquid part 

of the liquid- vapor film is solution of system : 

. . .dp ^ ^dp 

with A-^^^^=-7i+72Pu.o and A-^^^^ = -74 P,... • 

Quantity r is defined such that r = c// a/ Ap; , where l/r is a reference 
length and 73 = At. Solution of system (51) is 

P = P5 + Pi e"^^ + P2 e^^ (15) 

where boundary conditions at 2; = and h yield the values of pi and p2 : 

(72 + 73)pi + (72 - 73)P2 = 71 - 72P6 {S2) 
"^'"(73 - 74)pi + e'''^(73 + 74)P2 = -74Pb- 

The liquid density profile is a consequence of Eq. (15) when z € [0, ft.]. 
Taking Eq. (15) into account in Eq. (13) and gb{p) = {cf/2pi){p — pb)^ in 
linearized form for the liquid part of the interlayer, we get 

n(p.) = -^PiP2. (16) 

pi 

By identification of expressions (14), (16) and using (5*2), we get a relation 
between h and pb- We denote finally the disjoining pressure by Ii{h). 
Due to the fact that pb — pi,^"^ the disjoining pressure reduces to 

n{h) = — ^ [(71 - 72P;)(73 + 74)e''^ + (72 - 73)74P;] x 
Pi 

[(72 + lz)im - (71 ~ 72P;)(73 - 74)5"''^] 
[(72 + 73)(73 + 74)e''^ + (73 - 74)(72 - 73)e~'''"] 
Let us notice an important property of mixture of van der Waals fluid and 
perfect gas where the total pressure is the sum of partial pressures of com- 
ponents:^* at equilibrium, the partial pressure of the perfect gas is constant 
through the liquid-vapor-gas interlayer -where the perfect gas is dissolved 
in the liquid. The disjoining pressure of the mixture is the same than for 
a single van der Waals fluid and calculations and results are identical to 
those previously obtained. 
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5. Motions along a liquid nanolayer 

When the liquid layer thickness is small with respect to transverse dimen- 
sions of the wall, it is possible to simplify the Navier- Stokes equation which 
governs the flow of a classical viscous fluid in the approximation of lubrica- 
tion. -'^^ When h <^ L, where L is the wall transversal characteristic size, i) 
the velocity component along the wall is large with respcx't to the normal 
velocity component which can be neglected ; ii) the velocity vector varies 
mainly along the direction orthogonal to the wall and it is possible to neglect 
velocity derivatives with respect to coordinates along the wall compared to 
the normal derivative ; in) the pressure is constant in the direction normal 
to the wall. It is possible to neglect the inertial term when Re <^ L/h {Re 
is the Reynolds number of the flow). Equation of Navier-Stokes is not valid 
in a liquid nanolayer because the fluid is strongly inhomogeneous and the 
elastic stress tensor is not scalar. However, it is possible to adapt the ap- 
proximation of lubrication for viscous flows in a liquid nanolayer. We are in 
the case of long wave approximation: e = h/L <C 1. Wc; denote; the velocity 
by V = {u, V, w) where {u, v) are the tangential components. In the approxi- 
mation of lubrication we have : e = sup {\w/u\ ,\w/v\) <^1. The main parts 
of terms associated with second derivatives of liquid velocity components 
correspond to d^u/dz^ and d^v/dz^. The density is constant along each 
stream line (p = div'V = 0) and iso-density surfaces contain the tra- 
jectories. Then, du/dx, dv/dy and dw/dz have the same order of magnitude 
and e ~ e. As in Rocard model, we assume that the kinematic viscosity co- 
eflicient v = K2/ p depends only on the temperature.^*^ In motion equation, 
the viscosity term is (1/p) div ay = 2v [ div D + D grad {Ln(2«;2)} ] ; 
D grad{Ln {2k2)} is negligible with respect to div D. In both lubrication 
and long wave approximations the liquid nanolayer motion verifles 



r grad[ Ato(p) - A Ap] = z^AV with AV ~ 



dz'^ ' dz'^ ' 



(17) 



In approximation of lubrication, the inertial term is neglected and Eq. (17) 
separates into tangential and normal components to the solid wall. As in 
equilibrium, the normal component of Eq. (17) is 

d 

— [Ho{p)-XAp]=0 Ho{p) - XAp = ij.o{pb)- 

To each value pb (different of liquid bulk density value pi of the plane 

interface at equilibrium) is associated a liquid nanolayer thickness h. We 
can write ndp) — XAp = r}{h), where ry is such that -qih) = HoiPb)- For 
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one-dimcnsional motions colinear to the solid wall (direction io and velocity 
uio), the tangential component of Eq. (17) yields : 

grad Mo(P6) = lo ^ di,^=^d^- (18) 

A liquid can slip on a solid wall only at a molecular level. The sizes of 
solid walls are several orders of magnitude higher than slipping distances 
which are negligible and kinematic condition at solid walls is the adherence 
condition (z = u — Q). From the continuity of fluid tangential stresses 
through a liquid- vapor interface of molecular size and assuming that vapor 

viscosity stresses are negligible, we obtain [z = h, =^ — = 0). Conse- 

c) d f \. \ 
quently Eq. (18) implies vu= ^r— ^ ^r— ( - — hz\ . The mean spatial 

dpb dx \2 J 

velocity u of the liquid in the nanolaycr is u = u dz; previous com- 

_ /i2 ' _ ° _ . 

putations yield v u = — — grad Po{Pb) with u = m io. Let us remark 

that 

dpojpb) _ dpo_ dpb ^ ^ J_ dPjpb) dpb dh 
dx dpb dh dx pb dp dh dx 
The pressure Py^ in the vapor bulk is constant along flow motions and 



T-r/,x r./ X 1 duoipb) 1 dli(h) dh 

n(/i) = Py^ - P{pb); consequently, we get — = — — — and 

dx Pb dh dx 



X6 u= y gradn(/i). (19) 
where Xb = Pb^ is the liquid kinetic viscosity. 

Eq. (19) yields the mean spatial velocity of the isothermal liquid nanolayer 
as a function of the disjoining pressure gradient. Like as the disjoining 
pressure depends on the nanolaycr thickness, the mean flow velocity is a 
function of thickness variations along the flow. 
Taking into account that in the liquid nanolayer p ~ p^, then 



/ p dz\ u:^ pu 

,Jo j Jo 



dz 



and the mean spatial velocity corresponds also to the mean velocity with 

respect to the mass density. 

In shallow water approximation, the equation of continuity yields 

= 
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and we obtain Eq. (14) in ref. (^^) associated with /i— perturbations : 

dh , _ „ 
— + h divu = 0. 
at 

Thanks to Eq. (19) an equation for /i— perturbations is : 

Eq. (20) is an equation of diffusion in parabolic structure with a good sign 

of diffusion coefficient associated with stabihty when '^^[^'^ < 0.^^ 

an 
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